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WHAT IS CLAIMED IS: 

1 . A method for vascular elastography comprising: 

providing pre-tissue-motion and post-tissue-motion images in 
5 digital form of a vessel delimited by a vascular wall; said pre-tissue-motion 
and post-tissue-motion images being representative of first and second 
time-delayed configuration of said vessel; 

partitioning at least portions of both said pre-tissue-motion 
and post-tissue-motion images into corresponding data windows; 
10 approximating a trajectory between said pre-tissue-motion 

and post-tissue-motion images for corresponding data windows; and 

using the trajectory for each data window to compute a strain 
tensor in each data window. 

15 2. A method as recited in claim 1 , further comprising using 

said strain tensor in each data window to create an elastogram of at least 
part of said vessel. 

3. A method as recited in claim 1 , wherein said pre-tissue- 
20 motion and post-tissue-motion images are radio-frequency (RF) images. 

4. A method as recited in claim 3, wherein said pre-tissue- 
motion and post-tissue-motion images are part of a sequence of radio- 
frequency (RF) images. 

25 

5. A method as recited in claim 1 , wherein said pre-tissue- 
motion and post-tissue-motion images are issued from magnetic 
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resonance imaging (MRI), optical coherence tomography (OCT), 
brightness mode (B-mode)or Doppler-based ultrasound modality imaging. 

6. A method as recited in claim 1, wherein providing pre- 
5 tissue-motion and post-tissue-motion images in digital form of a vessel 

includes inducing tissue compression or dilatation on said vessel. 

7. A method as recited in claim 6, wherein inducing tissue 
dilatation on said vessel is achieved by cardiac pulsation. 

10 

8. A method as recited in claim 1 , wherein said strain tensor 
is the full strain tensor in at least one of said data windows. 

9. A method as recited in claim 8, wherein said full strain 
1 5 tensor is computed from three-dimensional or two-dimensional ultrasound 

data. 

10. A method as recited in claim 8, further comprising using 
said full strain tensor to compute the Von Mises (VM) coefficient in each 

20 data window. 

1 1 . A method as recited in claim 1 , wherein approximating a 
trajectory for each said data window includes using a Lagrangian speckle 
model estimator (LSME). 

25 
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12. A method as recited in claim 1, wherein a trajectory is 
approximated in each said data window using zero-order and first-order 
terms of a Taylor-series expansion, yielding: 
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[Tr] is a translation vector 

[LT] is a linear geometrical transformation of coordinates 
(x f y, z) represents the new position of a point (x 0 , yo, Zo). 



10 



15 



13. A method as recited in claim 12, wherein using said 
trajectory for each said data window to compute a strain tensor in each 
data window includes performing a non-linear minimization for each data 
window Wjj by computing a transformation [LT] providing the best match 
between each Wy of said pre-tissue motion image and a corresponding 
window Wjj in said post-tissue motion image. 

14. A method as recited in claim 12, further comprising 
computing the full strain tensor s having the components sy wherein: 

ey(t) = ^MO+Aj^t)] 
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(u, v, w) being displacement vector in the Cartesian coordinate system. 



15. A method as recited in claim 14, further comprising 
5 determining an elastogram providing a distribution of each component of 
the deformation matrix A and of the strain tensor e. 



16. A method as recited in claim 15, further comprising 
computing a pressure gradient between said pre-tissue-motion and post- 

10 tissue-motion images; said pressure gradient being used in determining 
said elastogram. 

17. A method as recited in claim 15, further comprising 
computing the Von Mises (VM) £ coefficient in at least some of said data 

15 windows as: 

4 = {|[( e xx - e yy) 2 +( £ yy "S 2Z ) 2 +(8 ZZ -S xx ) 2 +6(8^+8^+8^) 
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18. A method as recited in claim 17, further comprising 
determining a composite elastogram providing a distribution of the VM 
coefficient in at least some of said data windows. 

5 19. A method as recited in claim 17, further comprising: 

providing pressure gradient a resulting from blood flow 
pulsation of said vessel when said pre-tissue motion and post-tissue 
motion images are taken; and 

computing the elastic modulus in at least some of said data 

10 windows as: 




20. A method as recited in claim 12, wherein using said 
trajectory for each said data window to compute a strain tensor in each 
15 data window includes solving the following minimization equation: 

MIN |l(x(t 0 ) ? y(t 0 ) 5 z(to))-I Lag (x(t 0 + At),y(t 0 + At),z(t 0 +At))f 

where = [Tr;LT(:)] for data window Wy for augmented vector (;) and 
matrix vectorisation (:) 

I Lag (x(t 0 + At),y(t 0 + At),z(t 0 +At)) is the Lagrangian speckle image 
20 (LSI) defined as the post-tissue motion image 
l(x(t 0+ At),y(t 0 -hAt) ? z(t 0+ At)) numerica n y compensated for tissue 
motion. 
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21. A method as recited in claim 20, wherein solving said 
minimization equation includes using a minimization algorithm. 

5 22. A method as recited in claim 21, wherein said 

minimization algorithm is the regularized nonlinear minimization 
Levenberg-Marquardt (L&M) minimization algorithm. 

23. A method as recited in claim 12, wherein using said 
10 trajectory for each said data window to compute a strain tensor in each 
data window includes solving in a region of interest represented in both 
said pre-tissue-motion and post-tissue-motion images characterized by p x 
q pixels: 

Vi ! Xl yi Vi ! x, :*■ Vi V* Vi \ 

Ix 2 *2 r x 2 y2 !x 2 z 2 l x 2 ■■' l z 2 *2 Iz 2 Y2 Iz 2 Z 2 \ 2 



_* x pxq X P*q ^ x pxq ^Pxq ^pxq ^pxq ' ' ' ^Zpxq X p*q ^Zpxq ^pxq \y«q Z p«) I Zpcq 

15 for each corresponding said pixels in said pre-tissue motion and post- 
tissue motion images in digital form; 
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where t x =x(0,0,0,t); t y = y(0,0At); t z =z(0 5 0,0,t); and 
It = (l Lag (x(t + dt), y(t + dt), z(t + dt)) - l(x(t), y(t), z(t))) . 

24. A method as recited in claim 1, wherein providing pre- 
5 tissue motion and post-tissue motion images in digital form includes 

collecting longitudinal and cross-sectional radio-frequency (RF) data from 
said vessel. 

25. A method recited in claim 1 for endovascular 
10 elastography (EVE). 

26. A method as recited in claim 24, wherein providing pre- 
tissue-motion and post-tissue-motion images includes acquiring 
intravascular RF images using a catheter. 

15 

27. A method as recited in claim 26, wherein acquiring 
intravascular RF images using a catheter includes sequentially sweeping 
an ultrasound beam over a predetermined angle. 

20 28. A method as recited in claim 1 for non-invasive vascular 

elastography (NIVE). 

29. A method as recited in claim 25 for non-invasive 
microvascular elastography (MicroNIVE). 

25 

30. The use of the method from claim 1 for predicting risks 
of vascular tissue rupture or vascular aneurysms. 
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31 .The use of the method from claim 1 for phenotyping in 
animal models using genetic or cloning technologies. 

5 32. The use as recited in claim 31 wherein said model is 

hypertension (HT). 

33. The use of the method from claim 1 for in vivo 
measurements. 

10 

34. A system for vascular elastography comprising: 

an ultrasound system for acquiring pre-tissue motion and 
post-tissue motion radio-frequency (RF) images of a vessel; said pre- 
tissue motion and post-tissue motion images being representative of first 
15 and second time-delayed configuration of said vessel; 

a controller, coupled to said ultrasound system, i) for 
receiving said pre-tissue motion and post-tissue motion RF images, ii) for 
digitizing said pre-tissue motion and post-tissue motion RF images, iii) for 
partitioning both said pre-tissue motion and post-tissue motion RF images 
20 into corresponding data windows, iv) for approximating a trajectory for 
each said data windows; and v) for using said trajectory for each said data 
window to compute a strain tensor in each data window; and 

an output device coupled to said controller to output 
information related to said strain tensor in each data window. 



25 
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35. A system as recited in claim 34, wherein said controller 
further includes an analog-to-digital acquisition board for digitizing said 
pre-tissue motion and post-tissue motion images. 

5 36. A system as recited in claim 34, wherein said ultrasound 

system includes an ultrasound instrument, coupled to said analog-to- 
digital acquisition board. 

37. A system as recited in claim 36, wherein said ultrasound 
10 instrument includes a scanhead. 

38. A system as recited in claim 37, wherein said scanhead 
includes an array ultrasound transducer. 

15 39. A system as recited in claim 37, wherein said scanhead 

includes a single-element oscillating transducer. 

40. A system as recited in claim 36, wherein said ultrasound 
instrument includes a catheter having a tip and a transducer provided at 
20 said tip. 

41 .A system as recited in claim 36, wherein said ultrasound 
instrument is in the form of an ultrasound biomicroscope for non-invasive 
microvascular elastography (MicroNIVE) measurement. 



25 
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42. A system as recited in claim 36, wherein said ultrasound 
instrument is coupled to said analog-to-digital acquisition board via a 
radio-frequency (RF) pre-amplifier. 

5 43. The use of the system recited in claim 34 for predicting 

risks of vascular tissue rupture or vascular aneurysms. 

44. The use of the system from claim 32 for phenotyping in 
animal models using genetic or cloning technologies. 

10 

45. The use of the system from claim 44 wherein said model 
is hypertension (HT). 

46. A system for vascular elastography comprising: 

15 means for providing pre-tissue motion and post-tissue 

motion images in digital form of a vessel; said pre-tissue motion and post- 
tissue motion images being representative of first and second time- 
delayed configuration of said vascular vessel; 

means for partitioning both said pre-tissue motion and post- 
20 tissue motion images into corresponding data windows; 

means for approximating a trajectory for each said data 

windows; and 

means for computing a strain tensor in each data window 
using said trajectory for each said data window. 



25 



